The Astrophysical Journal, in press 

Preprint typeset using I4T]^X style emulateapj v. 02/07/07 



THE MAGNETAR XTE J1810-197: VARIATIONS IN TORQUE, RADIO FLUX DENSITY AND PULSE PROFILE 

MORPHOLOGY 

F. CAMILO, 1 I. COGNARD, 2 S. M. RANSOM, 3 J. P. HALPERN, 1 J. REYNOLDS, 4 N. ZIMMERMAN, 1 E. V. GOTTHELF, 1 D. J. HELFAND, 1 

P. Demorest, 5 G. Theureau, 2 and D. C. Backer 5 

Received 2006 October 22; accepted 2007 March 1 9 



o 
o 



OS 



(N 
> 

in 
oo 
\o • 
o ■ 

o ■ 

i ' 
o ■ 

a: 



ABSTRACT 

We report on 9 months of observations of the radio-emitting anomalous X-ray pulsar XTE J 1 8 1 0-1 97 starting 
in 2006 May using the Nancay, Parkes, GBT, and VLA telescopes mainly at a frequency of 1 .4 GHz. The torque 
experienced by the neutron star during this period, as inferred from a measurement of its rotational frequency 
derivative, decreased by 60%, although not in a steady manner. We have also observed very large ongoing 
fluctuations in flux density and pulse shape. Superimposed on these, a general diminution of flux density and 
a broadening of the pulse profile components occurred nearly contemporaneously with a decrease in torque 
of about 10% that took place in late 2006 July over an interval of 2 weeks. After a slight increase in average 
flux density, since 2006 October the flux density has continued to decline and the pulse profiles, while still 
varying, appear more uniform. In addition, a simultaneous observation of the pulsar with the Chandra X-ray 
Observatory and the GBT allows us to show how the X-ray and radio profiles are aligned. We discuss briefly 
the implications of these results for the magnetospheric currents in this remarkable object. 
Subject headings: pulsars: individual (XTE J1810— 197) — stars: neutron 



1. INTRODUCTION 

The transient anomalous X-ray pulsar (AXP) XTE J1810- 
197 (P = 5.54s) was discovered in ear ly 2003 when its X- 
ray luminosity increased ~ 100-fold (Ibr ahim et alJ 120041) 
compared to the quiescen t state maintained for > 24 years 
dHalpern & Gotthel i 120051) . Initial X -ray observations re- 
vealed unsteady spin-down with P« 10" 11 (implying a sur- 
face magnetic dipole field strength B s» 2 x 10 14 G) that var- 
ied by up to a factor of 2 within a few months of the out- 
burst. Such rot ational instability has been observed in other 
magnetars (see I Wood s & Thompson 20061 for a review of the 
properties of magnetars). 

The origin of a radio sou rce positionally coincident with the 
pulsar (Hal pern et alJl2005l) was clarified by the detection of 
strong, narrow, highly linearly pola rized radio pulses once per 
stellar rotation ( Camilo et al. 2006). This emission has several 
unique characteristics. First, it is at least one order of magni- 
tude brighter than in the late 1990s, when it was undetected; 
presumably the radio source turned on following the X-ray 
outburst. Just as the X-ray flux has no w decayed to a level 
appro ximating the historical low state ( Gotthel f & Halpera 
120071) . the radio emission can be expected to c ease "shortly". 
In the initial observations. [Camilo et al.l d2006l) also noted that 
the radio flux density of XTE J18 10-197 varied on time scales 
of approximately a day in a manner inconsistent with inter- 
stellar scintillation, implying large day-to-day radio luminos- 
ity changes not seen in ordinary pulsars. The "average" pulse 
profiles were also found to change in a manner seemingly 
inconsistent with variations observed in other neutron stars. 
Lastly, the magnetar has a very fiat radio spectrum, and has 
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been detected at higher radio frequencies than any other pul- 
sar. 

These unusual characteristics of XTE J18 10— 197 presum- 
ably reflect different physical conditions in its magnetosphere 
compared both with those of ordinary radio pulsars and of 
persistent magnetars, which so far h ave shown no evide nce of 
magnetospheric radio activity (e.g., Burgay et al. 2006|). The 
study of this radio emission thus provides a new window into 
the coronae of magnetars. We summarize here some results 
from our monitoring program of XTE J1810-197, focusing 
on rotational, flux density, and pulse profile evolution. 

2. DATA ACQUISITION, ANALYSIS, AND RESULTS 
2.1. Observations 
2.1.1. Nancay 

We have observed XTE J1810-197 at Nancay since 2006 
June 1 (MJD 53887) with typical integration times of 15- 
60min, using th e Berkeley-Orleans-Nancay (BON) coher- 
ent dedispersor (ICognard & Theureau! 120061) based on a 
Serendip V spectrometer 6 . The Nangay radio telescope has 
a gain of 1.5 KJy" 1 and a system temperature of 47 K at 
1.4GHz in the direction of XTE J1810-197, which this 
meridian-type telescope can track for 1 hr each day. We moni- 
tored the pulsar nearly every day until 2006 October, and have 
done so every 2-3 days on average since then. The last obser- 
vations reported on here are from 2007 January 27. Dedis- 
persion of a 64 MHz band centered on 1398 MHz is done co- 
herently into sixteen 4 MHz channels using a 64-node com- 
puter cluster, with each data stream then folded every 2 min 
at the pulsar's predicted period. In Figure Q] we show the 
daily folded pulse profiles obtained in the linear horizontal 
polarization (with dipole orientation parallel to the ground) af- 
ter excising significant levels of radio frequency interference 
(RFI). Regional telephone relay and radio transmitters use 
the vertical polarization, leading to much more severe RFI in 
that channel. However, t he radio emission f rom XTE J1810- 
197 is highly polarized (ICamilo et alJl2007h and fortuitously 

6 See |http://seti.berkeley.edu/casper/projects /SERENDIP5/ 
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FIG. 1 . — Daily pulse profiles from XTE J1810— 197 recorded in the lineal' 
horizontal polarizati on at N ancay at a frequency of 1.4 GHz across a band- 
width of 64 MHz (§ 12.1. It . each displayed with 2048 bins (the first is from 
2006 June 1 and the last from 2007 January 27). We aligned these profiles us- 
ing an ephemeris c ontain ing 1 1 frequency derivatives that resulted in feature- 
less residuals (see § 12. 2\ . Each profile is labeled with the effective integration 
time (i.e., after the removal of RFI) in minutes to its left and the observing 
date (MJD) to its right. In some cases minor artifacts remain, due to incom- 
pletely excised RFI. A high-resolution version of this figure can be obtained 
at http://www.astm. Columbia. edu/^fernando/XTEJ 1810-197Nancay.eps. 

is much stronger in the direction parallel to the ground at 
Nancay, so that the vertical signal is in any case usually very 
small (the parallactic angle change during the maximum inte- 
gration is a negligible ±2.°5). 

Since 2006 July we have used a second observing system in 
parallel with BON, where square-law detectors sampled every 
2 ms record the total power for each polarization channel in a 
50 MHz band centered on 1410 MHz. These data are flux- 
calibrated by the firing of a calibrated pulsed noise diode at 
the beginning of each observation. However, this detector's 
sampling rate is not stable enough for pulsar timing purposes, 
for which we use the accurately time-tagged BON records. 

We have used data from both systems to obtain flux den- 
sity estimates for the daily profiles. The square-law detected 
outputs suffer less from RFI (because of their higher cen- 
tral frequency and narrower bandwidth) and are accurately 
flux-calibrated. On the other hand, removal of RFI in these 
records is complicated by the lack of frequency resolution. 
For BON data we first carefully removed RFI in the time- 
frequency plane and then used the radiometer equation, along 
with the known telescope noise characteristics, to convert the 



rms noise fluctuation in the off-pulse part of the profiles to 
a Jansky scale. We determined this conversion factor for the 
horizontal polarization and applied it to both channels. We 
then added the corresponding profiles from both polarizations 
and, by integrating the area under each profile, finally mea- 
sured the pulsed flux density. On most days the flux densities 
estimated from the two systems agree to within ~ 25%, and 
hereafter we use only the BON values. 

It is evident from Figure Q] that the pulse profiles are not 
stable, particularly for the 2 months following late 2006 July 
(MJD 53935 is July 19). This presents obvious difficulties 
for obtaining pulse times-of-arrival (TOAs) by following the 
usual prescription of cross-correlating a "standard profile" 
with each day's profile. Instead, we obtain TOAs from the 
maximum of a parabola fitted to one pulse profile component 
that, we assume, regardless of variations, corresponds to the 
same longitudinal fiducial point on the pulsar. We choose for 
this fiducial point the left-most of the (typically) two compo- 
nents that make up the main cluster of emission. This choice 
is motivated by the facts that, within our large and densely 
sampled set of observations, (1) the longitudinal separation 
between the corresponding peaks (0.08P) remains essentially 
fixed even as the profiles vary enormously (e.g., see the pro- 
files for MJDs 53929 and 53977 in Fig.Q] for which the fidu- 
cial points are effectively the peaks of, respectively, the high- 
est and lowest clearly visible pulse components), and (2) the 
angular separation between the left-most peak and the "pre- 
cursor" visible on some days 0.3 in pulse phase preceding it 
(which is generally much fainter than the main components, 
but not always: see, e.g., the profiles for MJDs 53925 and 
53951) also remains constant within the uncertainties. Also 
commending this method is the empirical consideration that 
it works well, with only a handful of TOAs rejected because 
their pre-fit residuals are particularly large, owing presum- 
ably to "awkward" pulse shapes. The TOAs thus obtained 
for timing purposes (§ 12.21 ) also match (up to a constant offset 
of w 25 ms) those obtained via the standard cross-correlation 
method for data acquired before mid 2006 July (when profile 
variation was relatively small). 

2.1.2. Parkes 

The first detection of radio pulsations from XTE J18 10-197 
took place at the ATNF Parkes telescop e on 2006 March 17 , 
with confirmation following on April 25 dCamilo et al.l l2006). 
Since then we have observed the pulsar there at frequencies 
ranging from 0.7 to 8.4 GHz using a variety of spectrometers. 
Unless otherwise noted, in this paper we use only 1.4 GHz 
data, obtained with the analog filterbanks as described by 
ICamiloet all (120061) . 

The 1.4 GHz Parkes data set is very sparse by compari- 
son with Nancay's, but is particularly useful here for two 
purposes: to extend the overall timing span by more than 1 
month, and to perform consistency checks on the timing so- 
lution using two independent sets of TOAs. Because the vast 
majority of the Parkes data used here were obtained before 
late 2006 J uly (wh en pulse shapes were relatively stable; see 
Fig. [T]and § 12.1.11 ). the standard cross-correlation method that 
we used to obtain these TOAs works well enough. 

2.1.3. GBT 

We have observed XTE J1810-197 at the NRAO Green 
Bank Telescope (GBT) since 2006 May 2 at frequencies span- 
ning 0.3^-2 GHz. Very few of these data are at 1.4 GHz. 
Starting in 2007 January, after the monitoring rate decreased 
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FIG. 2. — Changing pulse shapes of XTE J1810— 197 from GBT observa- 
tions. All profiles were aligned with the polynomial used for Fig. [T] The 
labels list the observation date (MJD ), spectrometer us ed, and effective in- 
tegration time in minutes. Spigot I Kaplan et al. 2005) data were recorded 
across a 600 MHz band centered at 1 ,95 GHz, while a 134 MHz band cen- 
tered on 1.85 GHz was used for BCPM lBacker lit"aTlll997l) . Note the low- 
level emission near phases 0.8 on MJDs 53857, 0.6 on 53955, 0.65 on 53957, 
and 0.95 on 53963. See also Fig.|5] 

substantially at Nancay (§ 12.1.11 ). we increased GBT mon- 
itoring at 1.9 GHz, with data collected as described by 
ICamilo eTa l. (2006). The TOAs derived from these data, 
using the standard cross-correlation method with pulse pro- 
files that vary little, are used to complement Nancay TOAs. 
Also, using the large gain, bandwidth, and frequency agility 
of GBT, we have obtained very high-quality pulse profiles at 
a number of widely-separated frequencies that illustrate the 
remarkable time-variation in pulse shapes displayed by this 
magnetar. Some examples of this collection at a frequency of 
1.9 GHz are shown in Figure [2] 

2.1.4. VIA 

We have used a sub-array of the NRAO Very Large Ar- 
ray (VLA) to measure with high precision the flux density of 
XTE J1810-197 at 1.4GHz on 13 days between 2006 Febru- 
ary 28 and September 5. During a typical observation we ac- 
cumulate 40 min of on-source visibilities in continuum mode, 
with a bandwidth of 100 MHz. Each observation begins with 
a w 5 min scan of a flux calibrator, either 3C286 or 3C48. 
The remainder of the observation consists of three « 2 min 
scans of the phase calibrator 181 1-209 interspersed with two 
« 20 min on-source scans. The recorded visibilities are cal- 
ibrated and imaged using standard procedures in the Astro- 



nomical Image Processing System. Due to variations in the 
array configuration, the number of antennas used, and the 
amount of observing time allocated, the errors in our flux den- 
sity measurements range between 0.2 mJy and 1.0 mJy. 

2.1.5. Chandra 

XTE J1810-197 was observed with the Chandra X- 
ray Observatory on 2006 September 10-11, using the 
back-illuminated S3 chip of the ACIS-S CCD detector in 
TIMED/VFAINT mode with a sub-array readout that pro- 
vided a time resolution of 0.441 s, and deadtime of 9%. Data 
were acquired continuously over an interval of 8.4 hr start- 
ing at 19:40 UT with no background contamination. All pho- 
ton arrival times were corrected to the solar system barycen- 
ter (TDB) using the pulsar position given in Table Q] the 
JPL DE200 ephemeris and the axbary task in the Chandra 
Interactive Analysis of Observations. Photons falling within 
an aperture of radius 3" centered on XTE J1810-197 were 
extracted from the standard processing event Level 2 data 
files and a total of 7014 counts were accumulated in the 0.5- 
3.0keV energy range, chosen to maximize the source signal- 
to-noise ratio. 

The barycentered photon arrival times were folded at epoch 
MJD 53989.0 (midnight TDB on September 11) into 10 bins 
using P = 5.540362s determined from a simultaneous radio 
observation at the GBT (see § 12.2. 11 1. 

2.2. Timing 

We have used the available TOAs (§§ l2.1.ll42.1.3l l. along 
with the TEMPO 7 software, to obtain the timing solution 
that describes the rotational history of XTE J18 10— 197 dur- 
ing 2006 May-2007 January. Because of the ch angin g pulse 
shapes and heterogeneity of the available data (§ 12. Il l, partic- 
ular care must be taken with this process. 

First we used Nancay TOAs to obtain a solution beginning 
in 2006 June. Starting with a handful of TOAs and an ini- 
tial solution fitting only for rotational phase and frequency 
v = l/P, we subsequently added v as needed. We then added 
one TOA at a time, paying particular attention to its pre-fit 
residual and corresponding pulse shape, as well as the post- 
fit residuals and fit parameters. In this way we eliminated 
seven TOAs, from a set of 144, that had large (> 100 ms) pre- 
fit residuals; these could usually be ascribed to particularly 
unusual pulse shapes. After approximately 1 month it is no 
longer possible to obtain featureless timing residuals without 
fitting for higher frequency derivatives. After this we added 
the Parkes TOAs, fitting also for an offset between both TOA 
sets to account for an arbitrary alignment between the respec- 
tive fiducial points. The magnitude of this offset (21 ± 15 ms; 
unless otherwise stated, the uncertainties used in this paper 
all represent 1 a confidence levels) is as expected, consider- 
ing the differ ent metho ds used to obtain the respective TOAs 
(see §§ 12.1.11 and 12.1.2b . Finally, when measuring v in piece- 
wise fashion for 2007 January, we complemented the Nancay 
TOAs with those from the GBT (§|2T3). 

In general it is difficult to estimate TOA uncertainties re- 
liably for our data, owing to changing pulse shapes and the 
principal method used to obtain TOAs (see § 12.1. 11 1. For the 
early Parkes data these are typically in the 5-10 ms range. The 
uncertainty of early Nancay TOAs is more often 10-20ms. 
Following the general weakening of the pulsar, the broadening 
of its profile components, and the huge variability displayed 

7 See http://www.atnf.csiro.au/research/pulsar/tempo/ 
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FIG. 3. — Timing residuals for XTE J1810— 197. Black dots co rresp ond to 
Nan9ay TOAs, while open circles represent Parkes TOAs (see § 12. It . Top: 
Residuals for a model that contains only the rotation phase, frequency v and 
v (see Tableff), showing a clear cubic residual trend. Bottom: Residuals for a 
model with an additional polynomial term, v = 9.40 ±0.06 X 10~ 21 s~ 3 . Note 
the change in vertical scale by a factor of about 10, and the remaining quartic 
trend. 



by the pulse shapes (Fig. [TJ, we estimate that the TOA uncer- 
tainties are commonly around 25-50 ms but can reach values 
approaching 100 ms on occasion. As noted above we have ex- 
cluded from the timing fits a few TOAs with such large pre-fit 
residuals (and estimated uncertainties). Because of these dif- 
ficulties, we assigned a uniform weight to the TOAs and, as 
a result, we do not obtain a \ 2 figure of merit for the timing 
fits. Nevertheless, the TOA uncertainty estimates noted here 
are in keeping with the corresponding fit rms residuals. 

We obtained a phase-connected timing solution that encom- 
passes all data beginning with the confirmation TOAs from 
April 25 (we comment below on use of the discovery TOAs 
from 2006 March 17). The residuals from a quadratic fit to 
pulse phase are shown in the top panel of Figure [3] and the re- 
spective solution is given in TableQ] An independent solution 
using only Parkes TOAs results in similar parameters. The cu- 
bic trend of large amplitude visible in the residuals points to 
very significant unmodeled behavior, expected of young pul- 
sars in general (w hose rota tion is "noisy") and of XTE J 1810- 
197 in particular (Ibrahim et al. 2004). While v clearly varies 
with time, the bottom panel of Figure[3]shows that it does not 
do so smoothly. A total of 1 1 frequency derivatives are re- 
quired to "whiten" the residuals, which then have an rms of 
18 ms. 

In order to extract from this record of unsteady rotation a 
quantitative measure of the varying torque acting on the neu- 
tron star, we have measured v as a function of time. We have 
done this in a piece-wise fashion, performing quadratic fits to 
phase for TOA spans long enough that the nominal resulting 
fractional uncertainty in v was no more than 3%, but short 
enough that no trends were seen in the residuals. In practice 
this implied individual segments of very nearly 30 days each. 
So as to better sample the variation in i), we did this with 
Nancay data (supplemented in 2007 January by GBT TOAs) 
while stepping through the TOAs in offsets of 15 days. A fit 
to all the good Nancay TOAs using 1 1 frequency derivatives 
also yields v(t) (without associated uncertainties) consistent 
with the discrete measurements. 

Owing to the lack of Nancay observations and the sparser 
data from Parkes and GBT, this was not straightforward be- 



TABLE 1 

Timing parameters for XTE J 1 8 10—1 97 



Parameter Value 

Right Ascension 18 h 09 m 51?087 

Declination -19°43'5i;'93 

DM(cnT 3 pc) 178.0 

Epoch (MJD) 54000.0 

v(s-') a 0.180493891(6) 

v (10~ 13 s~ 2 ) a -2.53(1) 

MJD range 53850-54127 

NOTE. — The celestial coordinates were held 
fixed at the values ob tained from VLBA observations 
iHelfand et al. 2007), and the DM was held fixed at 
the value obtained from simultaneous 0.7 and 2.9 GHz 
observations (Camilo et al. 2006). 
a These two parameters are sufficient to obtain 
a phase-connected solution encompassing the MJD 
range, but do not fully describe the rotation of the 
neutron star. They are non- statio nary and, strictly, not 
predictive. See Fig.[5]and jj |2.2l for more details. 

fore 2006 June 1. We opted for performing timing fits us- 
ing multi-frequency data from both Parkes and GBT, thereby 
gaining some additional leverage in our solutions. At this 
time, the radio emission at most frequenc ies was usually dom - 
inated by a single narrow component (see lCamilo et aU 2006), 
so that TOA artifacts were relatively small. In this manner, we 
obtained a v measurement for late April-late May. Finally, we 
also used the 2006 March 17 Parkes discovery observations in 
a timing fit. We appear to have maintained phase connection 
with the confirmation observations 39 days later (the pre-fit 
residual between the discovery TOAs and the solution com- 
puted with the first month of data following confirmation is 
only 0.16P). However, we have established with Nancay data 
that cubic trends typically become significant for > 30 days, 
and we have therefore obtained this earliest measurement of 
v also by computing the difference in barycentric frequen- 
cies between 2006 March 17 and April 25. The two v values 
differ by only 2.7 a of the (much smaller) uncertainty of the 
phase-connected fit, and this agreement allows us to extend 
our record of v by 1 additional month. 

The run of v measured over 9 months is shown in the top 
panel of Figure |4] The overall trend is one of increasing v (as 
could be inferred from the top panel of Fig. [3) but the way 
in which this occurs is far from steady. In the 3 months prior 
to mid 2006 July, v varied in a relatively steady and grad- 
ual fashion, from -3.3 x 10" 13 s" 2 to -3.0 x 10" 13 s" 2 . Then, 
in late July, it changed to -2.7 x 10~ 13 s" 2 over a span of 15 
days — the implied change in torque is -2 x 10 32 dyn cm (us- 
ing a stellar moment of inertia of 10 45 g cm 2 ), which is larger 
than the torque powering three-quarters of all known ordi- 
nary pulsars! Regardless of the detailed mechanisms that pro- 
duce radio pulses in XTE J1810-197, it should perhaps not 
be surprising if other observational properties of the magnetar 
should have changed around this time. 

In fact, as the middle panel of Figure [4] shows, the peak 
flux density of XTE J18 10-197 has dramatically decreased 
since about that time. Interestingly, while the period-averaged 
flux density has also decreased compared to its average value 
before mid July (bottom panel of Fig. HI, until about 2006 
October it did so by a smaller factor and continued to fluctu- 
ate greatly from day to day (when average flux densities are 
available from at least two of Nancay, Parkes or VLA within 
a day of each other, they are consistent within expectations 
given the inherent variations). These two observations can be 
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FIG. 4. — Frequency derivative and flux densities at 1.4 GHz for 
XTE J1810— 197. Top: Frequency derivative obtained from a variety of tim- 
ing fits each spanning about 30 days. Solid error bars denote fits obtained 
with Nancay data alone (excepting the last one, which also contains some 
GBT data), while dotted error bars correspond to fits that contain Parkes and 
in some cases also GBT data. For fits containing Nancay data the error bars 
represent nominal 2a confidence levels; for all others, 4<r. In late 2006 July 
(indicated by the dashed line at MJD 53940 in the lower panels), v increased 
at a large rate, the flux densities re ache d a local minimum, and the pulse pro- 
files changed in character. See § 12.21 for details. Middle: Daily pea k flux 
density of the profile component used to obtain Nancay TOAs (see § 12.1.11 
and Fig.fT). Bottom: Daily period-averaged flux density from Nancay (black 
dots) and Parkes {open circles), and continuum VLA flux density {triangles 
with error bars). 

understood by inspection of the profiles shown in Figure [T] 
from late July to mid September (MJD — 53994), the daily 
profiles of XTE J1810-197 tended to be composed of two (or 
more) significant peaks, each much broader than the one peak 
generally prominent before mid July (with typical full-width 
at half-maximum « 0.04P, versus about half that value be- 
forehand), and with greater pulse-shape variance than before. 
This may also explain why the Nancay timing residuals for 
each 1 -month fit after mid 2006 July are about 24 ms rms, 
twice as large as the typical corresponding value before then. 

Since 2006 October, the Nancay profiles appear to have 
varied less, and to be composed mainly of one broad peak 
(Fig. OX although on some days the trailing peak is still rec- 
ognizable (as it is more often in higher quality GBT data). 
After another relatively large increase over the month of 
September, v has continued a steady increase at approxi- 
mately the average rate for the 9-month span of our obser- 
vations, 7.5 x 10~ 16 s" 2 per day. 

2.2.1. X-ray and radio pulse alignment 

On 2006 September 10-11 we observed XTE J1810-197 
for 7.2 hr with the GBT at 1.9 GHz, start ing 1 hr after the 
beginning of the Chandra observation (§ 12.1.51 ). We have 
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FIG. 5.— Bottom: X-ray pulse profile (0.5-3.0keV) of XTE J1810-197 
from an observation with Chandra on 2006 September 10-11, folded with 
barycentric P = 5.54 0362 s as determined from a simultaneous radio obser- 
vation (see § 12. 1.51 . Phase zero on this plot corresponds to MJD 53989.0 
(TDB). Top: Radio profiles from a simultaneous GBT observation, starting on 
September 10, 20:50 UT. After 6 hr (bottom trace), the pulse profile abruptly 
changed to a different configuration for 30min (middle trace), after which 
it changed yet again (top trace). The relative integrated flux densities of the 
three profiles are 1.4, 1.0, 2.9, from bottom to top. The middle peak of the 
bottom radio profile arrives a t phas e 0.17 on the plot, and all four profiles 
are absolutely aligned (see 8 12.2. It . Note also the alignment between the 
left-most peak of the top profile and the first notch in the bottom profile. 



used the contemporaneous TOAs thus obtained to measure the 
phase offset between the X-ray and radio pulses. The folded 
X-ray profile is shown in Figure [5] with the first phase bin 
chosen to begin at midnight on September 1 1 (TDB). After 
translating the TOAs for the first radio profile (see Fig. [5]) con- 
verted to infinite frequency to the solar system barycenter, a 
TEMPO fit yields a phase offset for the fiducial point of the 
radio profile (here its middle peak) of 0.167 ±0.006 in the 
Figure. The uncertainty includes a component due to the fit 
and a slightly smaller contribution from the uncertainty in the 
dispersion measure (DM) of XTE J1810-197. 

Within the larger uncert ainty imposed by the Chandra reso- 
lution (0.08P; see § I2.1.51 l and the relatively small number of 
pulsed X-ray counts, the main component of the radio profile 
on this day arrives at the same time as the peak of the X-ray 
pulse (Fig. [5}. 

Interestingly, the radio profile changed in both shape and 
flux during the Chandra observation — twice, within a span 
of 30min (Pig. [5j. This does not modify the conclusion in the 
above paragraph, but helps answer a question concerning the 
radio flux and profile variations observed in XTE J18 10— 197: 
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they can occur suddenly (observed at a resolution of ~ 10 s). 
We have "caught in the act" at least nine such changes at 
Nancay, GBT and Parkes, all since mid 2006 July. As this 
corresponds to some 150 hours of observing time, we can es- 
timate that such transitions occur on average every ~ 15hr, 
at the present epoch. No X-ray bursts or significant changes 
in the X-ray flux or pulse shape were seen at the times of the 
radio transitions indicated in Figure [5] The Chandra count 
rates during the three radio pulse states — 0.235 ± 0.003, 
0.213 ±0.011, and 0.218 ± 0.007 s" 1 , respectively — are un- 
corrected with the large changes in radio flux density and 
pulse shape. 

3. DISCUSSION 

All young radio pulsars experience rotational instabilities, 
observed as a continuous quasi-stochastic wandering of the 
pulse phase ("timing noise") or as discontinuous spin-up 
events ("glitches"). These are thought ultimately to be driven 
by the unsteady transfer of angular momentum from the inte- 
rior superfiuid to the crust of the neutron star. XTE J1810- 
197 displays large amounts of timing noise-like behavior; de- 
spite very large changes observed in v on short time scales, all 
observed rotational changes have been continuous to within 
the available resolutio n and we have not observe d any clear 
glitch-like behavior. [Arzoumani an et ail (Q994) quantified 
timing noise as the time-magnitude of its cumulative con- 
tribution over a time interval t to the cubic term in a Tay- 
lor series expansion of rotational phase, i.e., vt^ /(6v). For 
XTE J18 10— 197 this amounts to 120s over 277 days, a huge 
amount that far surpasses anything observed in radio pul- 
sars and is also greater than an extrapolation based on the 
notion that the magnitude of timing noise is proportional to 
P dArz oumania n et al]|1994l) . However, this level of timing 
noise is not unprecedented for magnetars. 

Six AXPs, including XTE J1810-197, and two soft-gamma 
repeaters (SGRs) — i.e., the majority of all known mag- 
netars — have been observed with phase connection main- 
tained for at least a few months (see iKaspi et all 1200 It 
Gotthelfetai] l2002t iGavriil & Kaspil l2002t IWoods et all 
2002t IKaspi & Gavriill 12003b iGavriil & Kasol 12004 The 
SGRs 1806-20 and 1900+14 usually have the largest amounts 
of timing noise, roughly one order of magnitude above the 
level we observ e in XTE Jl 8 10-1 97, with v varying by fac- 
tors of up to 4 (IWoods et a l. 2002). The rotational stability of 
AXPs is generally g reater (although at least two have been ob - 
served to glitch; see lKaspi et alj2003tlKaspi & Ga vriil 2003), 
and four are substantially quieter than XTE J1810— 197. The 
remaining one, IE 1048.1-5 937, was relatively quiet dur- 
ing 2005-2006 dKaspil 120071) . but previously displayed far 
greater rotational insta bility than observed in XTE J18 10-197 
(IGavriil & Kaspil 120041) . Therefore, while noisier than most 
AXPs, XTE J1810-197 may not be particularly remarkable 
in this respect. What is unusual is that we can now track its 
rotation via nearly daily observations, in greater detail than is 
possible for other AXPs. 

Following its discovery with the Rossi X-ray Timing 
Explorer, a phase-connected solution was obtained for 
XTE J1810-197 spanning 2003 January-September, dur- 
ing which v varied between -6.7 x 10~ 13 s" 2 , averaged over 
8 months, and -3 .8 x 10~ 13 s" 2 over the last 2 months 
(llbrahimetal.112004 . The very large initial magnitude may 
have reflected in part transients associated with a putative 
glitch at the time of the X-ray outburst in early 2003. In the 
succeeding 2.5 years, until the radio detection, phase connec- 



tion was no longer possible owing to the sparse sampling, but 
individual period measurements 6 months apart established, 
e.g., that the average v was -1 .7 x 10~ 13 s~ 2 be tween late 2003 
and early 2004 (see iGotthelf & Halpernll2005l) . The value of 
v = -1 .4 x 10~ 13 s" 2 in early 2007 is the largest ever measured 
for this source, for which the historical range can be summa- 
rized as v = (-2.6 ± 1.2) x 10~ 13 s" 2 spanning late 2003-early 
2007. Importantly, although v has increased almost monoton- 
ically during our radio observations (Fig. 0), it cannot have 
been increasing during the entirety of the past 3 years (see also 
IGotthelf &^ Halpern 2007). Periods of decreasing torque such 
as we have observed recently must have been interspersed 
with at least one epoch of substantially increasing torque. If 
the torque for XTE J1810-197 were to continue decreasing at 
the rate in 2007 January, by mid 2007 the pulsar would stop 
spinning down altogether; presumably this rate will soon de- 
crease. 

In ordinary pulsars, the observed rotational instabilities are 
mainly presumed to be caused by angular momentum trans- 
fer processes internal to the neutron star. In the magne- 
tar model, the external energy stored in the magnetosphere, 
partly released through reconfiguration and decay of the mag- 
netic field and currents, may also lead to erratic variations 
in the torque acting upon the star, for example via a flux 
of Alfven waves and particles resulting from magnetically- 
induced (sudden large -scale or persistent sma ll-scale) crustal 
seismic activity (e.g., Thom pson et alj|2000l). It is unclear 
whether any such particular models (see also iDuncanlfeOOll) 
can explain the magnitude and time scale of the variations in 
torque currently observed in XTE J 1 8 1 0-1 97, 4 years after the 
X-ray outburst. In any case, substantial variations in magne- 
tospheric plasma densities and/or currents would likely have 
implications for other observable properties of the pulsar. 

An interesting comparison is with PSR B 193 1+24, an ordi- 
nary radio pulsar during intervals of 5-10 days, with P = 0.8 s 
and Z? « 2 x 10 12 G, which abruptly shuts off for 25-35 days i n 
a pattern that repeats quasi-periodically (Kr amer et alJl200a) . 
During these turned-off p eriods, the torque is only two-thirds 
of its turned-on value. iRramer et alj ([2006) conclude from 
this extraordinary behavior that the occasional presence of 
plasma leads to radio emission and its flow provides the ex- 
tra braking torque. The plasma density calculated from the 
torque differen ce is 100 esu cm" 3 , in agr eement with the co- 
rotation value dGoldreich & Julian! fl969l) . In XTE J1810- 
197 the radio emission has not shut off, but it did dimin- 
ish and change markedly in character somewhat abruptly in 
late 2006 July, apparently coincident with a huge reduction 
in torque (i> changed from -2.96 x 10" 13 s" 2 on MJD 53933 
to -2.65 x 10" 13 s" 2 on MJ D 53949; see §|Z2lan d Fig.gji. A 
calculation such as that by Kram er et alJ ((2006) would sug- 
gest for XTE J1810-197 an "extra" plasma density of about 
300 esu cm" 3 prior to the large torque decrease, compared to 
the source's Goldreich-Julian density of 1350 esu cm" 3 . Im- 
plicit in this calculation is the notion that the radio emission 
originates from open field lines, which has not been proven for 
XTE J18 10-197. For a magnetar, there is additional plasma 
due to non-axisymme tric large-scale magnetospheric currents 
(iThompson et al.l2002l) . In any case, the pulse profile and flux 
variations observed in XTE J1810-197 occur on time scales 
that are apparently too sho rt to be explained by changes in th e 
closed field lines (see, e.g.. lBeloborodov & Thomp son 2007). 
Also, interpretation of polarimetric data for XTE J1810- 
197 appears consis tent with emission from open field lines 
(Cami lo et afll2007l) . 
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If indeed the measured changes in torque are caused by 
variations in the magnetospheric plasma density in loca- 
tions relevant for the production of coherent radio emis- 
sion, then it should be no surprise that the pulse profiles of 
XTE J1810— 197 change so remarkably: even the "small" 
torque changes reflected in Figure|4]are enormous by the stan- 
dards of ordinary pulsars. And certainly, the profile variations 
in XTE J1810-197 (Figs. [T] E] and EJl are observationally dis- 
tinct from the "mode changing" of some ordinary pulsars, 
where the average pulse profile suddenly chan ges between 
two o f a small set of different configurations (e.g., Batt el et all 
1982); XTE J1810— 197 displays a much greater variety. How- 
ever, the root cause of mode changing is not well understood 
and, as we have established that XTE J18 10— 197 profiles can 
change suddenly (Fig. |5}, some link may exist between these 
two phenomena. Other behavior also appears extraordinary, 
such as the broader pulse components typically observed af- 
ter MJD 53940 even as the phase separation between com- 
ponents remains fixed (Fig. [TJ and § |2.2t and the polarim etric 
properties remain largely unchanged (Cam ilo et al .1 12007). 

Compared to such extreme radio variability, the X-ray spec- 
trum and pulse profil es of XTE Jl 8 10—197 chan ge slowly as 
the flux decays (e.g., Gotthelf & Halpern 2007b, on a much 
longer time scale (~ 1 yr) than the fluctuations in torque re- 
ported here. The alignment of the peaks of the radio and X- 
ray pulses suggests that the footpoints of the active magnetic 
field lines on which radio emission is generated are also the 
locations of concentrated crustal heating that is responsible 
for the enhanced X-ray emission, at least at the higher ener- 
gies that apparently come from a relatively small area (see 
Gotthelf & Halpern 2007). However, even if the energetic 
particle bombardment that heats the surface hot spot fluctu- 
ates on time scales of less than 1 day, the absence of correlated 
X-ray variability on similar time scales indicates that most of 
the X-ray luminosity originated in deeper crustal heating at 
the time of the X-r ay turn on, or from mor e gradual decay of 
the magnetic field dEichler& Chenglll989t) . 

Among persistent AXPs it is not altogether clear what is 
the relationship between observed radiative properties (such 
as pulse shapes and fluxes) and rotational evolution (such 
as variations in torque). Pulse profiles can change follow- 



ing glitches (e.g. .lKaspi et al.ll2003l). but otherwise appear to 
be fairly stable ( Gavriil & Kaspi 2002). Quite apart from 
the pheno menon of X-ray b ursts (for bursts in XTE J1810- 
197, see IWoods etalJ 12005). some AXPs have shown sub- 
stantial variability in X-ray flux (e.g.. lGavriil & K aspi 2004). 
Nevertheless, with the p ossible exception of IE 1048.1-5937 
dGavrii l & Kaspi 2004), there has been no reported correla- 
tion between X-ray flux and spin-down rate . For XTE J 1 8 1 0— 
197, on the other hand, there is an observed correlation be- 
tween radio flux density and torque over a period of 9 months 
(Fig.|H, although we know neither whether this correlation is 
causal nor whether it holds over longer periods of time. 

With observations of XTE Jl 8 10-197, we now find our- 
selves in the curious position of relying upon energetically 
insignificant radio pulsations exhibiting a remarkably diverse 
phenomenology to illuminate rather more energetic events on 
the magnetar. Attempting to understand key aspects of the ra- 
dio observations may lead to a deeper understanding of both 
magnetars and radio pulsar emission. 
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